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FOREWORD 


This  rcporc  was  prepared  as  part  of  a research  program  at  the  Naval 
Weapons  Center  on  the  subject  of  Ignition  of  Solid  Propellant  Rockets^ 
sponsored  under  Naval  Ordnance  Systems  Comaand  Task  Assignment  ORD  331- 
00I/200-1/UR024-02-02.  Previous  reports  include  an  extensive  sumnary  of 
the  then  ^current  status  of  various  theories  of  solid  propellant  ignition 
(NAVWEPS  Report  8987,  NOTS  TP  3954),  and  extended  analyses  of  solid  phase 
theory  (KVC  TP  4618)  and  heterogeneous  theory  (NWC  TP  4864).  The  present 
report,  supported  in  addition  by  Ballistic  Research  Laboratories  MIPR 
RN-516-75,  is  the  third  in  a series  of  three  reports  on  a Unified  Theory 
of  Solid  Propellant  Ignition.  It  presents  some  results  of  a parameter 
study  obtained  through  use  of  the  theory  and  computer  program  developed 
in  the  first  two  parts. 


This  report  has  been  prepared  for  timely  presentation  of  information. 
Because  of  the  continuing  nature  of  research  in  this  area,  refinements 
may  be  made  in  the  future. 
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(U)  /I  Unified  Theoru  of  Solid  Propellant  Ignition. 
Par:  Computer  Solutions  by  H.  H.  Bradley*  Jr. 

China  Lake*  Calif.*  Naval  Weapons  Center*  December 
N.  1975;.  28  pp.  (NWC  TP  5618*  Part  3*  publication 

TSO^iVIFIED). 

4^  Computer  solutions  are  obtained  for  several 
ignition  models  involving  either  surface  or  gaseous 
exothermic  chemistry.  Surface  regression  and  con- 
vection of  products  from  the  surface  are  found  to  be 
limiting  for  the  surface  reaction  of  gaseous  oxidizer 
and  solid  fuel.  Fast  gas  phase  reactions  cannot  in 
general  he  distinguished  from  surface  and  solid  phase 
reactions  by  the  experimental  determination  of  the 
dependence  of  ignition  tine  on  igniting  flux;  a more 
usable  technique  nay  be  the  effect  of  concentration 
on  ignition  regimes.  Slov  gas  phase  reactions  would 
be  best  studied  by  photographic  observation  of  the 
developing  gaseous  reaction  zone.  Further  analyses, 
coupled  with  judiciously  chosen  numerical  solutions* 
offer  a reasonable  path  to  the  solution  of  the  re- 
maining unresolved  questions. 
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lOTRODUCTION 


Ignition  of  combustible  solid  materials  depe*)ds  upon  many  interrelated 
phenomena,  including  chemical  kinetics  and  heat  and  mass  transport.  A 
general  ignition  theory  would  d*^al  with  all  processes  simultaneously  but 
is  too  complicated  for  tractable  solution  or  presentation  of  results  of 
calculations.  When  mass  diff us ional  processes  are  neglected,  the  ensuing 
theory  is  called  the  thermal  ignition  theory.  Heat  losses  to  the  gas 
phase  may  be  included  but  are  usually  of  lesser  importance.  Several 
reviews  (e.g. , Price,*  and  Merzhanov  and  Averson^)  have  been  published  ‘ 
dealing  with  thermal  ignition  theory  and  with  various  mathematical  tech* 
niques  of  solving  the  equations  either  numerically  or  by  appropriate 
approximate  methods.  The  approximate  solutions  agree  exceptionally  well 
with  the  more  exact  numerical  solutions  and  encompass  a variety  of  models 
and  boundary  conditions,  so  that  further  exploration  in  this  field  does 
not  seem  fruitful.  Less  restricted  models  (e.g.,  Williams,^  and  Linan 
and  Crespo**),  including  surface  reactions  and  diffusional  processes,  have 
been  presented.  The  broader  field,  in  which  diffusion  and  convection 
play  an  important  role,  has  not  yet  been  so  exhaustively  studied.  The 
Importance  of  the  various  processes  appears  in  proper  perspective  if  it 
is  recognized  that 

1.  Kinetics  and  thermal  transport  govern  time  to  incipient  Ignition 
(temperature  runaway). 

2.  Diffusional  mass  transport  establishes  limits  for  incipient 
ignition. 

3.  Convection  controls  the  transition  from  incipient  ignition  to 
steady-state  combustion. 

The  success  of  approximate  analytical  methods  of  solutions  in  the 
case  of  thermal  Ignition  models  rests  entirely  with  the  great  sensitivity 
of  the  chemical  reaction  rate  to  temperature  over  the  range  of  activation 


* Price,  E.  W. , and  others.  "Theory  of  Ignition  of  Solid  Propellants,** 
AIAA  J,  Vol.  4,  No.  7 (July  1966),  pp.  1153-81. 

^ Merzhanov,  A-  C.,  and  A-  E.  Averson.  **The  Present  State  of  the 
Thermal  Ignition  Theory;  An  Invited  Review,**  COMBUST  AND  FLAME,  Vol.  16 
(1971),  pp.  89-124. 

’ Williams,  F.  A.  ''Theory  of  Propellant  Ignition  by  Heterogeneous 
Reaction,**  AIAA  J,  Vol.  4,  No.  8 (August  1966).  pp.  1354-57. 

**  Linan,  A.,  and  A.  Crespo.  **An  Asymptotic  Analysis  of  Radiant  and 
Hypergolic  Heterogeneous  Ignition  of  Solid  Propellants,**  COMBUST  SCI  AND 
TEaiNOL,  Vol.  6 (December  1972),  pp.  223-32. 
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energy  encountered,  thereby  localizing  the  chcolstry  and  permitting  divi- 
sion of  the  ignition  process  into  an  initial  inert  heating  phase  followed 
by  a cheoical  phase  of  ntpldly  increasing  intensity.  With  decreasing 
activation  energy,  the  chemistry  becomes  more  diffuse  in  time  and  position, 
leading  to  a reduction  in  accuracy  and  ultimate  breakdown  of  approxima- 
tions. In  systems  containing  large  diffusional  or  convectional  effects, 
a situation  similar  to  low  activation  energy  may  prevail.  Indeed,  it  has 
been  noted'  that  interpretation  of  experimental  results  to  yield  activa- 
tion energies  nay  be  in  error  when  mass  transport  effects  are  ignored. 

While  the  characteristics  of  Ignition  of  systems  not  containing  im- 
portant diffusional  or  convectional  effects  arc  quite  well  known,  the 
relative  paucity  of  information  regarding  the  other  end  of  the  spectrum 
leads  to  complications  in  the  analysis  of  experimental  results  and  makes 
difficult  the  choice  between  the  various  ignition  mechanisms.  Two  methods 
suggest  themselves  for  approaching  the  problem.  The  first  involves  ex- 
tension of  the  techniques  of  approximation  used  in  thermal  ignition 
theory,  especially  that  of  dividing  the  process  Into  stages,  each  of 
which  may  be  treated  separately.  The  second  would  be  the  development  of 
an  entirely  different  concept  of  approximation  tailored  to  the  convection- 
diffusion  problem.  As  in  the  development  of  thermal  ignition  theory, 
accurate  numerical  solutions  are  required  against  which  to  test  approxi- 
mations. 

This  report  is  the  third  in  a series  of  which  the  first*  presented  a 
unified  ignition  model  containing  as  many  of  the  diverse  ignition  situa- 
tions as  feasible;  the  second^  described  a computer  program  for  solving 
the  equations  developed  in  Reference  6.  The  current  report  presents  some 
preliminary  results  of  computer  calculations.  It  is  not  the  purpose  of 
this  report  to  provide  a handbook  of  solutions  spanning  a wide  range  of 
parameters.  An  exhaustive  or  even  comprehensive  parameter  study  using 
numerical  methods  would  be  prohibitive  in  terras  of  time  and  money;  rea- 
sonably accurate  approximate  analytical  methods  are  needed  for  such  pur- 
poses. Future  reports  in  the  series  will  be  addressed  to  the  study  of 
such  analyses. 

The  framework  for  the  present  study  was  developed  in  Reference  6, 
which  contains  complete  details.  To  place  the  current  effort  in  proper 
perspective,  the  principal  assumptions  by  which  the  equations  of  the 
ignition  i.iode’i  were  derived  from  conservation  laws  are  listed  without 
comment: 


' Rosner,  D.  E.  "Convective  Diffusion  as  an  Intruder  in  Kinetic 
Studies  of  Surface  Catalyzed  Reactions,"  AIAA  J,  Vol.  2,  No.  4 (April 
1964),  pp.  593-610. 

^ Naval  Weapons^  Center.  A Unified  Theory  of  Solid  Propellant  Ignition^ 
Part  i.  Development  of  Mathematical  Models  by  H.  H.  Bradley,  Jr.  China 
Lake.  Calif.,  N^^C,  August  1974.  36  pp.  (NWC  TP  5618,  Part  1.) 

’ Naval  Weapons  Center.  A Unified  Theory  of  Solid  Propellant  Ignition^ 
Part  2.  Computer  Program^  by  H.  H.  Bradley,  Jr.  China  Lake,  Calif.,  NWC, 
August  1974.  28  pp.  (NWC  TP  5618,  Part  2.) 
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1*  One  space  dimension 

2.  No  cross-coupling  of  forces  and  fluxes 

3.  No  radiation  feedback  from  gas  phase  reaction  zone 

4.  No  momentum  transport  (constant  pressure) 

5.  Temperature-independent  properties 

6.  Semi-infinite  gas  and  solid  phases 
Validity  of  Chapman-nubesin  approximation 

8.  No  subsurface  gasification  in  solid 

9.  Adequacy  of  global  kinetics 

To  simplify  computer  programming,  some  additional  approximations 
were  used  in  Reference  7,  namely; 

10.  Constant  heat  flux  boundary  condition,  hence  no  absorption  of 
incoming  stimulus  by  evolved  gases 

11.  Chemistry  far  removed  from  equilibrium,  so  that  rate  effects  are 
governing 

12.  Sensible  heat  small  compared  to  heats  of  reaction  and  phase 
change 

13.  Adequacy  of  the  binary  diffusion  approximation 

14.  Equal  molecular  weights  of  all  species 

In  order  to  restrict  the  amount  of  computational  work  still  further, 
the  general  model  was  specialized  to  four  simple  models  to  be  described 
In  subsequent  sections.  Also,  some  additional  limitations  were  imposed 
on  the  calculations: 

15.  Completely  opaque  solid  phase — no  subsurface  absorption  of  in- 
coming radiant  flux 

16.  No  solid  phase  reactions 

17.  Uniform  initial  temperature — equal  in  the  solid  and  gas 

Assumption  16  was  made  because  of  the  availability  of  comprehensive 
studies  for  the  case  of  solid  phase  reactions.®"^® 


SELECTION  OF  CASES  FOR  STUDY 


Ignitablc  systems  may  be  placed  into  two  categories  according  to 
whether  all  the  ingredients  participating  in  the  ignition  chemistry 


* Linan,  A.,  and  F.  A.  Williams.  "Radiant  Ignition  of  a Reactive 
Solid  with  In-Depth  Absorption,"  COMBUST  AND  FLAME,  Vol.  18  (February 
1972),  pp.  85-97. 

® Bradley,  H.  H. , Jr.,  "Theory  of  Ignition  of  a Reactive  Solid  by 
Constant  Energy  Flux,"  COMBUST  SCI  A1  1 TECHNOL,  Vol.  2 (August  1970), 

pp.  11-20. 

Linan,  A.,  and  F.  A.  Williams.  "Theory  of  Ignition  of  a Reactive 
Solid  by  Constant  Energy  Flux,"  COMBUST  SCI  AND  TECHNOL,  Vol.  3 (April 
1971).  pp.  91-98. 
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arise  frara  the  propellant  Itself  or  whether  at  least  one  reactant  is  pro- 
vided as  an  elem  .it  of  the  igniter  in  addition  to  its  thermal  stimulus. 

A different  breakdovm  is  provided  by  the  site  (gas  or  surface)  of  the 
key  reaction  which  initiates  steady  state  combustion.  Solid  phase  bulk 
reactions  are  excluded  from  this  study  by  Assumption  16.  We  have  then 
the  possibility  of  four  subtheories  of  the  unified  ignition  theory: 

1.  Case  A:  Surface  chemistry  between  a solid  fraction  and  an 

ambient  reactive  gas. 

2-  Case  B;  Surface  chemistry  between  one  pyrolyzed  and  one  un- 
pyrolyzed  solid  fraction.  This  is  Case  A with  the  ambient 
reactive  gas  provid'd  by  pyrolysis  of  the  solid  itself. 

3.  Case  C:  Cas  phase  reaction  between  a pyrolyzed  solid  fraction 

and  an  ambient  reactive  gas. 

4.  Case  D:  Gas  phase  reaction  between  two  pyrolyzed  ingredients 

initially  present  in  the  solid. 

In  addition,  a brief  study,  accompanied  by  computer  calculations, 
was  made  on  endothermic  and  exothermic  surface  pyrolysis,  since  one  or 
both  processes  are  basic  to  several  of  the  four  main  subtheories. 

In  all  four  cases,  the  stimulus  may  be  provided  by  radiant  dux,  as 
in  arc  image  experiments,  or  by  an  initial  temperature  difference  be- 
tween phases,  as  in  shock  tube  experiments.  Combinations  of  one  or  more 
of  the  simplified  cases,  as  well  as  other  possibilities  are  roost  likely 
as  ignition  mechanisms;  nevertheless,  it  is  believed  that  thorough 
understanding  and  correlation  of  the  characteristics  of  the  simpler 
cases  are  basic  to  the  solution  of  the  overaxl  ignition  problem. 

Comprehensive  studies  of  Case  A with  both  temperature  difference  and 
radiant  heating  have  been  reported.**  This  report  will  concentrate  on 
only  one  aspect  of  Case  A,  namely,  the  effect  of  convective  mass  flux, 
engendered  by  surface  reaction,  on  the  character  of  the  solution  for  a 
representative  range  of  parameters. 

No  studies  of  Case  B are  currently  available,  although  results  on  a 
related  model  with  exothermic  surface  reaction  rather  than  distributed 
solid  reaction  have  been  reported.*^  The  distinction  between  Case  B 
and  Reference  12  is  that  in  the  latter,  no  mechanistic  modeling  of  gas 
phase  processes  is  attempted,  boundary  condi^io.is  being  based  on  a 
reasonable  application  of  steady  state  relationships.  A moderate  effort 
was  expended  on  Case  B to  provide  orientation  for  future  studies. 


**  uradley,  H.  H. , Jr.,  and  F.  A.  Williams.  "Theory  of  Radiant  and 
Hypergolic  Ignition  of  Solid  Propellants,"  COMBUST  SCI  AND  TECHNOL, 

Vol.  2 (August  1970),  pp.  41-52. 

*^  Baer,  A.  D. , and  N.  W.  Ryan.  "An  Approximate  but  Complete  Model 
for  the  Ignition  Response  of  Solid  Propellants,"  AIAA  J,  Vol.  6,  No.  5 
(May  1968),  pp.  872-77. 
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A presentation  of  Case  C with  radiant  stimulus  has  been  published. 
including  the  frequently  neglected  effects  of  convection  and  surface 
regression.  For  this  reason  Case  C will  not  be  covered  here;  computer 
calculations  were  limited  to  check  runs  to  establish  program  capability 
and  execution  times. 

Studies  of  Case  D are  seemingly  limited  to  instances  with  a tempera- 
ture difference  (relevant  to  shock  tube  work)  and  also  are  combined  with 
Case  C.***  Thus  they  apply  to  the  interesting  situation  in  which  two  oxi- 
dizers compete  for  the  fuel  vapors  in  gas  phase  reactions.  The  complexity 
of  Case  D.  owing  to  its  greater  number  of  parameters,  has  precluded  any 
comprehensive  study  in  this  report.  Some  representative  calculations 
have  been  made,  primarily  for  comparison  with  analogous  runs  of  Case  B. 


STUDIES  OF  INDIVIDUAL  CASES 


It  is  almost  universal  practice  to  use  dimensionless  parameters  in 
the  analysis  of  complicated  systems  of  equations.  The  required  number 
of  calculations  or  experiments  is  usually  significantly  reduced  if  the 
optimum  set  of  parameters  is  chosen.  An  optimum  set  is  one  from  which 
only  a few  meid>ers  have  an  important  influence  on  the  outcome  of  the 
calculations  or  experimental  results.  Often  such  an  optimum  set  cannot 
be  recognized  a priori,  especially  in  studies  Involving  new  or  little 
explored  fields.  In  such  instances,  a sparse  coverage  of  the  parameter 
space  renders  interpretation  of  results  in  terms  of  actual  physical 
quantities  almost  impossible.  Therefore,  the  choice  was  made  to  deal 
with  dimensional  equations  directly,  although  a dimensional  analysis 
was  carried  out  to  establish  the  minimum  number  of  dimensionless  param- 
eters required  to  describe  the  problem.  The  single  exception  was  Case  A. 
for  which  an  extensive  study  had  already  been  completed.^*  Otherwise, 
standard  values  were  assigned  to  those  physical  quantities  (Table  1) 
which  were  considered  constant  throughout  the  study,  while  a reference 
set  was  chosen  f ^ r the  other  parameters  about  which  to  conduct  a sensi- 
tivity study.  No  claim  is  made  for  the  relevance  of  the  reference  set 
to  a particular  physical  system;  the  choice  was  made  from  among  published 
values  and  from  other  studies. 


Kashivagi,  T.  **A  Radiative  Ignition  Model  of  a Solid  Fuel,” 

COMBUST  SCI  PSD  TECHKOL,  Vol.  8,  No.  5 and  6 (1974).  pp.  225-36. 

***  Kumar.  T.  K..  and  C.  E.  Hermance.  ‘‘Ignition  of  Homogeneous  Solid 
Propellants  under  Shock  Tube  Conditions:  Further  Theoretical  Development.” 

AIAA  J,  Vol.  9,  No.  8 (August  1971),  pp.  1615-20. 
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TABLE  !•  S^ndard  Physical  Propcrtlss  Us#d  Ui  i*^rc. 


Synbol 

Definition 

V4ltM 

Units 

Solid  doaslty 

• • • 

u^ 

B/c«’ 

Solid  specific  beat 

e 

on* 

. 3U0->) 

cal/g  •E 

Solid  cooductielty 

\ 

... 

laor^) 

csl/ca  *%  sec 

Solid  responsieicy 

Ape 

2.25#(10"*) 

cal/cn*  *K  sec^^^ 

Solid  penilctielty 

5.667(W>) 

g*/cn^  sec 

Cas  density  (initial) 

^20 

... 

1.016(10”*) 

g/cn* 

Gas  specific  heat 

c 

... 

3(10”‘) 

cal/g  *E 

Cas  conductivity 

^2 

... 

*O0”*> 

cal/ce  *K  sec 

Gas  mspoosivity 

r 

•2 

1.104(10”**> 

cal/cn’  •K  sec'^^ 

Cas  pemittiTity 

*2 

1.354(10-’) 

g*/cn^  sec 

Fressnre 

... 

... 

1 

atjs 

Cas  molecular  reigbt 

U 

... 

25 

Initial  ecsiperature 

^0 

... 

300 

•E 

Lewis  Bunber 

In 

*2 

1 

dimer  iionless 

PYROLYSIS 

Transient  pyrolysis  is  described  by  the  following  set  of  equations 
adapted  from  Reference  6. 

Field  equations: 

Solid  enthalpy 

5/},  a,  m db. 

1 s 1 

^10  10  1 

Gas  enthalpy 

3/)2  ^2  3^/)2  vn  ^^2 

Pyrolyzed  ingredient 

9y  ^ ^2  3^y  ”s  3y 

at  “ p|g  " P20  3*2 


(1) 


(2) 


(3) 
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Boundary  equations: 
Energy 


pio  " ho 


m Q 
r s 


(4) 


0 


0 


Pyrolyzed  Ingredient 


0 • 0 K - 


s s 


(5) 


where 


■ p,v  e*p(-E/jpr^) 

In  Che  absence  of  gas  phase  reactions  it  can  be  shown* ^ that  the  gas 


(6) 


phase  temperature  approaches  a spatially  constant  value  equal  to  the 
steady  state  surface  temperature  (T|)  and  given  by  solving  the  following 
steady  state  heat  balance  equations  for  the  surface: 


In  Reference  15 « a matched  asymptotic  expansion  vas  completed  on  the 
energy  equations  alone,  vith  the  result  that  the  surface  temperature 
was  found  to  approach  its  steady  state  value  asymptotically  from  below. 
Integration  of  the  current  equations^  including  raass  fiow  effects, 
reveals  a temperature  which  overshoots  the  steady  state  value  Tj  by 
2 to  20Zp  becoming  most  pronounced  for  exothermic  reactions  with  high 
activation  energy.  This  is  explainable  by  the  action  of  pyrolyzed  mass 
in  carrying  away  high- temperature  surface  material  in  increasing  amounts 
as  surface  regression  increases  toward  its  final  steady  state  value  as 
given  by  solution  of  Equation  7 and  8.  The  calculations  do  not  indicate 
any  effect  of  gaseous  mass  diffusivity  (manifest  in  Le)  on  the  shape  of 
the  curve  of  surface  temperature  vcrsus>  txce,  thus  implying  that  convec- 
tion, not  diffusion,  is  controlling  in  this  instance.  An  interesting 
and  unexplained  phenomenon  of  the  solution  of  Equations  1 through  6 is 
the  presence  of  small  amplitude,  low  frequency  damped  oscillations  in 
surface  tc-operature  for  endothermic  pyrolysis.  No  oscillations  were 
noted  for  cxothenaic  pyrolysis. 


Kindelan,  M. , an*  F.  \.  Williams.  "Theory  for  Endothermic  Gasifi- 
cation of  a Solid  by  a Constant  Energy  Flux,"  COMBUST  SCI  AND  TECHNOL, 
Vol.  10,  No.  1 and  2 (1975),  pp.  1-19. 


(7) 


(8) 
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CASE  A 


Tho  following  governing  equations  for  Case  A may  be 
by  reference  to  Part  1 of  this  report:^ 

Field  equations: 

Solid  enthalpy 

B/i 


dt 


a, 


m dhi 
s * 


"10 


B\|ii 


Gas  enthalpy 


Bt 


Ojo 

Oxidizer 

a 


m Bh^ 
s ^ 


^2  0 


■ Le 


2 


3<»> 


-s 


at 

Product 


PjO  ^20  ^"^2 


at 


Le 


-2 


Pro 

Boundary  conditions: 
Energy 
Bh 


^!!± 

P20  ^^^2 


1 1 
Pjo 


^20 


3 


Oxidizer 


n y 
s u 


0 

Product 


a r 

s 6 


Lc 


Lo 


P20 


ar, 

2 


2 , 

+ b, 

U|0  ■*  6 


^20* 


b-B- 

U 6 6 
^10 


where 


"s  “ “6  ■ exp(-2r^/RT^) 

Taking  constant  specific  heats  and  defining 


■ <’’i  - 


readily  vritten 

(9) 

(10) 

(ID 

(12) 

(13) 

(14' 

(15) 

(16) 
(17) 
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it  nay  be  shown  by  dimensional  analysis  that  the  solutions  of  Equations 
9 through  16  arc  expressible  in  the  functional  form 


V " ^1-  V 

(18) 

where 

(19) 

t - (</r;r5)t 

(20) 

(21) 

^ - C (a  /a  /p  )♦ 

(22) 

■ the  parameter  set: 

C,/Cj,  I».  Ej/HTj.  .ltd 

In  Reference  11,  advantage  Is  taken  of  an  Integral  solution  to  repre- 
sent the  surface  conditions  as  a function  of  time  with  no  regression  in 
terms  of  three  parameters: 

(23) 

6 - 

(24) 

(25) 

where  dimensionless  time  l is  defined  as 

T - Ky  + l)/Y"T.rVt  (26) 

, *6 

and  D replaces  the  expression  Pj 0^  * present  study 

of  Ciise  A,  the  same  parameters  are  used  for  convenience  in  comparing 
results.  The  additional  parameter  is  introduced  and  all  the  Equations 
9 through  16  are  solved  in  arriving  at  results.  Again,  only  results  for 
the  surface  are  presented,  otherwise  the  full  six  parameters  would  be 
required  to  describe  temperature  and  concentration  distributions.  The 
choice  i)^  * 1 represents  reaction  of  one  mass  unit  each  of  oxidizer  and 
fuel,  leading  to  evolution  of  two  units  of  product  and  to  a net  convec- 
tive flow  from  the  surface.  Assignment  of  smaller  values  of  would 

correspond  to  a more  severe  convective  effect  since  is  the  ratio  of 
oxidizer  mass  to  fuel  mass  reacted. 

Before  describing  the  results  of  the  parameter  study  it  is  convenient 
to  define  three  characteristic  temperatures  of  the  process.  First,  the 
steady  state  temperature  in  the  absence  of  external  heating  (g^  ■ 0)  and 
ignoring  convection  Is  shown  in  Reference  11  to  be  given  by 
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e - lAi 

s.^ 


(27) 


This  Is  the  dlncnslonlfss  tompcrnture  nt  which  the  rate  of  cheralcal  re- 
action of  oxidizer  just  balances  Its  rate  of  diffusion  to  the  surface. 
Second*  ve  have  the  temperature  at  whlch»  In  the  absence  of  diffusional 
or  convectTonal  effects*  the  applied  flux*  q * equals  the  chemically 
generated  surface  heat  ^ 

. In  y/(3  - In  Y)  (28) 

This  condition  obviously  cannot  exist  If  B < In  >.  Third*  there  Is  the 
steady  state  pyrolysis  tomper«iturc  Oj*  obtained  by  solving  the  pair  of 
Equations  7 and  8.  This  dimensionless  temperature  6|^  to  Indicate 

exothenilcity*  is  generally  greater  than  0^*  although  the  trend  rewrses 
at  higher  values  of  Oy.  U’e  also  note  here  that  the  steady  state  adl«i- 
batlc  ■ 0)  surface  temperature  with  regression  is  given  by  Equation  7 
because  at  steady  state,  as  noted  previously*  the  gas  pVase  temperature 
approaches  a space-uniform  value.  Hence  no  sensible  heat  is  lost  to  the 
gas  phase  and  all  the  surface  chemical  heating  goes  into  heating  the 
solid.  Under  these  conditions 


T 

ss  conv 


^0 


or  in  dimensionless  toms 


(29) 


0 conv  * 
ss 


(’•_ 


O /C  T 

W \ 0 


(30) 


Only  two  steady  state  solutions  exist  for  the  distribution  of  , namely* 
y^  » 0 and  \\  “ u^b^[/Ujn-  Corresponding  to  y,^  « 0 is  zero  concentration 
gradient,  mass  flux,  and  chemical  heating.  In  a number  of  numerical 
calculat ions,  external  heating  was  interrupted  (to  simulate  go/no-go 
experiments)  at  various  times  during  heating.  In  all  cases*  regardless 
of  the  state  of  the  system  at  the  time  of  interruption*  the  ensuing 
transient  led  ultinatuly  to  an  extinguished  state  in  which  y,^  and 
dl’i^/d^^  * 0 at  the  surface.  Tlie  go/no-go  criterion  of  ignition  is  thus 
seen  to  bo  inapplicable  to  ignition  corresponding  to  Case  A. 


Turning  now  to  Figure  1,  we  see  the  results  of  calculations  for  con- 
tinuous healing  over  a broad  range  of  values  of  parameters  a*  B*  and  Y- 
Numbered  points  on  the  figure  arc  keyed  to  run  numbers  of  Table  2*  which 
presents  other  quantitative  results.  Classical  ignition  curves  charac- 
terized by  a rapid  tonperature  rise  (Figure  2a)  occurred  only  in  Kuns  1 
and  2.  Even  those  exliibited  the  transition  to  zero  surface  chemical 
heating  rate  described  above  for  Interrupted  heating.  For  this  reason 
an  alternative  definition  of  ignition  time*  used  in  Reference  11  and 
represented  by  the  solid  curves  in  Figure  1*  was  employed.  This  was 
the  tine  of  occurrence  of  equal  applied  and  chemical  heating  rates* 
which,  without  convection,  was  often  close  to  the  time  of  rapid  tempera- 
ture rise  or  the  go/no-go  time.  With  the  single  criterion  chosen,  con- 
vective and  nonconvect ivc  times  were  nearly  equal*  but  the  range  of 
parameters  for  ignition  to  occur  was  reduced.  The  ircquality 

Gy(numerical)  < l/a  ^ (31) 
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/ 


NWC  TP  5618,  P#rt  3 


FIGURE  1.  Comparison  of  Ignition  Times  for  Case  A With  Surface 
Regression  (Points)  and  Without  Surface  Regression  (Solid  Lines). 
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TABLE  2.  Results  of  Case  Calculations. 


Run  no. 

a 

e 

o 

O 

0 ^ 
Y 

< 

1/a 

1 

0.1 

50 

10 

0.855 

0.86 

1.063 

10 

2 

1 

50 

10 

0.855 

0.86 

1.063 

1 

3 

5 

50 

3 

0.16 

0.17 

1.000 

0.2 

4 

10 

50 

1 

0.0482 

0.058 

1.000 

0.1 

5 

0.1 

30 

8 

1.591 

1.59 

1.680 

10 

6 

1 

30 

6 

0.855 

0.92 

1.126 

1 

7 

3 

30 

1 

0.083 

1.000 

0.2 

8 

0.1 

20 

6 

2.234 

a 

2.160 

10 

9 

1 

20 

2 

0.299 

0.32 

1.000 

1 

10 

5 

20 

1 

0.130 

0.32 

1.100 

0.2 

11 

0.1 

10 

2 

0.855 

0.86 

1.330 

10 

12 

1 

10 

1 

0.299 

0.35 

1.050 

1 

^ Calculated  by  Equation  28. 

From  numerical  solution. 

^ Calculated  from  Equation  7.8. 

Kot  reached. 

Is  still  required  for  ignition  as  in  the  nonconvect ive  case.  A new 
characteristic  associated  with  convection  is  the  presence  of  oscillations 
of  the  relaxation  t>T>e  (Figure  2b » super jjiposed  on  the  rising  portion  of 
the  curve  of  surface  temperature  versus  time  and  occurring  after  the 
temperature  has  risen  beyond  the  value  corresponding  to  0^.  Each  example 
of  oscill.itory  bchav'ior  is  associated  with  the  inequality 

< l/a  < (32) 

including  Runs  3,  4,  6.  7.  and  9.  Runs  1.  2.  5.  11,  and  12,  which  are 
near  the  periphery  of  the  parameter  space  in  Figure  1.  show  little  or  no 
tendency  toward  oscillation;  tempera tore- time  traces  resemble  the  lower 
portion  of  the  curve  of  Figure  2a.  No  ignition  occurs  in  Run  8 where  0^ 
is  never  reached  or  in  Run  10  for  which  Oy  > 1/a.  While  the  correlations 
arc  not  perfect,  they  suggest  that  the  interaction  between  chemical  reac- 
tion and  convection  alternately  consumes  the  gas  phase  reactant  at  the 
surface  and  blocks  it  from  entering  until  provided  by  the  slower  diffu- 
sional  process.  The  implication  of  this  convectional  blocking,  unique 
to  Case  A,  is  that  some  other  transport  mechanism,  not  included  in  the 
present  model,  must  bo  responsible  for  ignition  or  combustion  processes. 
The  nechanistn  is  probably  turbulent  mixing,  generated  in  real  situations 
by  cross  flow  of  gases  or  by  the  combustion  itself.  The  result  of  such 
nixing  would  be  to  reduce  the  effective  value  of  a,  thereby  extending 
the  range  of  other  parameters  for  which  ignition  occurs. 
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Characteristic  ignition  Ignition  with  oscillatory 

surface  tesperature 

FIGURE  2.  Typical  Dimensionless  Curves  of  Surface  Temperature 
History  During  Heating  by  Constant  Radiant  Flux  (Case  A). 


CASE  B 


There  are  two  rate  processes  in  Case  B:  (1)  endothermic  pyrolysis 

of  one  constituent  of  the  solid  propellant  to  produce  an  active  gaseous 
species,  and  (2)  exothermic  reaction  of  this  active  gas  in  a surface 
process  with  the  unpyrolyzed  fraction.  The  equations  solved  were: 

Field  equations: 


Enthalpy — same  as  Equations  9 and  10 


Oxidizer 

3*^2  3^2 

^20  °'*2  *^10  2 


Product 


10 


(33) 


(34) 
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Boundary  conditions: 
EnorgX 


Oxidizer 


^2  ^^2 


'20 


34', 


+ + a,  iQ 


u«’n  + 9, 


a 3r  y 

’ll  * *5^2  - ^ 5^7  347  117^  *2*5 


Product 


where 


^2  ‘'s 

0 ■ M y - Le  - — —, l>c™i- 

s 5 020  P‘5'2  PlO  5 t- 


m ^ + in. , 

s ^ 11 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


Wlien  the  sane  teopcracurc«  time,  and  distance  scales  are  used  (see 
Equations  17,  20,  21,  22),  dimensional  analysis  yields  the  following 
functional  solution: 


where 


^2^  “ ^1’  "2*  V 


the  parameter  set: 


(41) 


aj/a^.  c,/c^.  r^,  0^/c^T^,  E^/RT^, 

The  last  two  measure  the  relative  intensities  of  the  two  rate  processes 
and  the  stimulating  flux  at  the  initial  temperature  T^. 

Tlie  rapid  proliferation  of  parameters,  oven  in  the  nondimensional 
representation,  suggests  the  proliibitivc  number  of  c Iculations  which 
would  be  associated  with  a complete  parameter  study.  Even  the  more 
modest  approach  of  conducting  a sensitivity  study » which  ignores  inter- 
actions among  parameters,  requires  the  choice  of  a reference  set  of 
parameters  for  variation-  The  decision  was  to  make  a limited  study  of 
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the  effect  of  varying  the  physical  quantities  v^*  Vj ^ P. 

^ with  values  of  the  other  quantities  as  shown  in  Table  3 
The  variation  of  was  carried  out  in  two  ways:  (1)  independently, 

and  (2)  with  to  provide  a constant  reference  value  of  the  exothemic 
reaction  rate  between  fuel  and  pyrolyzed  oxidizer.  Likewise  E^^  was 
varied  independently  and  with  to  yield  a constant  reference  value  of 
oxidizer  pyrolysis  rate.  The  radiation  flux  level  was  varied  from  3 to 
30  cal/cm^sec,  pressure  effect  was  evaluated  from  0.1  to  10  atmospheres, 
and  Levis  number  was  varied  from  0.1  to  10.  The  effect  of  and  C|| 

was  found  by  keeping  their  algebraic  sun  equal  to  2000  cal/g,  thereby 
allowing  both  endothermic  and  exothermic  pyrolysis. 


TABLE  3.  Parameters  Used  in  Case  B Calculations. 


Parameter 

Reference  value 

Range 

30  kcal/mole 

25  - 35 

'*5 

7.26(10*’) 

1.35(10*^)  - 3.91(10^^) 

30  kcal/mole 

25  - 35 

''ll 

2.91(10"*) 

5.40(10**)  - 1.56(10*^) 

7 

10  cal/cin^sec 

3-30 

r 

0 

S 

2150  cal/g 

1850  - 2000 

^11 

-150  cal/g 

0 - 150 

Lc 

1 

0.1  - 10 

P 

1 atm 

0.1  - 10 

Three  chatacteristics  of  the  solution  were  examined:  (1)  the  time 

(t^)  to  attain  zero  surface  chemical  heating  (endothermic  heat  of  pyroly- 
sis equal  to  exothermic  heat  of  reaction),  (2)  time  (t2)  at  which  total 
surface  chemical  heating  equals  Incident  flux,  and  (3)  time  (tj)  of  occur- 
rence of  the  inflection  point  in  the  curve  of  surface  temperature  versus 
time  when  such  a ooint  exists.  For  the  parameter  values  examined,  t2 
and  agree  within  5%. 

The  effect  of  radiant  flux  level  on  and  t2  is  displayed  in  Fig- 
ure 3 for  the  reference  set  of  physical  parameters.  The  relationship  is 
similar  to  results  obtained  for  solid  phase  bulk  reactions  and  for  pure 
exothermic  surface  reaction.  An  essential  difference  from  Case  A is  the 
absence  of  convective  blocking  of  reactants  by  products,  although  the 
loss  of  energy  fr6m  the  surface  would  have  a marked  effect  during 
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r-yiTTim 1 — 1 rmiT 

100 

1 \ E 
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3 10  30 

q,  cal/cm^sec 

FIGURE  3.  Effect  of  Flux  Level  on  Ignition 
Times  tj  and  t2  for  the  Reference  Set  of 
Parameters  (Case  B). 

transition  to  steady  state.  No  quantitative  test  of  this  phenomenon  has 
been  made,  although  either  suppression  of  surface  regression  or  inter- 
rupted stimulus  runs  would  be  important. 

The  influence  of  propellant  kinetic  parameters  is  shown  in  Figure  A. 
The  effect  of  varying  either  activation  energy  with  its  pre-exponential 
factor  is  nearly  undetectable  in  terms  of  characteristics  (1)  and  (2); 
however,  there  is  no  inflection  point  unless  at  least  one  activation 
energy  'ixc?-.eds  an  approximate  value  of  30  kcal/mole.  When  either  activa- 
tion energy  is  varied  independently,  the  results  are  much  more  dramatic, 
ranging  from  very  short  ignition  times  to  no  ignition  at  all.  When  the 
pyrolysis  activation  energy  is  lower  than  the  reaction  activation  energy, 
no  ignition  of  type  (1)  or  (2)  was  observed  to  occur.  Physically,  this 
is  consistent  with  the  temperature-limiting  effect  of  the  surface  endo- 
therm,  which  prevents  the  attainment  of  a temperature  high  enough  to 
promote  the  exothermic  reaction  within  practical  time  limits.  Exposure 
of  the  solid  fuel  to  constant  temperature  reactive  pyrolysis  products 
might  lead  to  a situation  similar  to  an  "adiabatic"  ignition  were  it  not 
for  the  effect  of  convection  in  removing  any  developing  thermal  zone  from 
the  surface.  When  pyrolysis  activation  energy  is  higher,  the  potential 
for  reaction  occurs  prior  to  the  availability  of  both  species;  as  the 
pyrolysis  temperature  is  approached,  exothermic  reaction  occurs  rapidly 
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at  the  onset  of  pyrolysis.  Ignition  tiaes  within  the  ignition  regimes 
increase  with  either  activation  energy  when  such  energy  is  varied  inde- 
pendently of  its  associated  pre-exponential  factor. 


Eg.  E,,  KCAUMOLE 

FIGURE  4.  Effect  of  Activation  Energy  on  Ignition 
Time  t2  for  the  Reference  Set  of  Parameters  (Case  B) . 


The  effect  of  pressure  is  similar  to  the  effect  of  transport  proper- 
ties in  imposing  an  ignition  limit  but  not  otherwise  affecting  ignition 
tine  significantly-  At  a pressure  of  1/10  atmosphere,  characteristic  (2) 
was  not  established  so  that  ignition  would  not  occ^r.  At  10  atmospheres 
the  ignition  time  was  decreased  approximately  52.  The  result  is  antici- 
pated since,  by  the  Chapman-Rubesin  approximation.  is  constant,  and 
diffusional  processes  with  their  induced  convection  and  surface  cooling 
are  more  pronounced  at  low  pressures. 

Ignition  occurred  over  the  entire  range  of  Lewis  numbers  investigated. 
For  Le  « 0-1,  corresponding  to  low  mass-diffusion  rates,  the  ignition 
time  was  reduced  from  27  to  24  millisei *>nds,  while  for  Lo  * 10  an  increase 
to  28  milliseconds  was  noted  bas^c  again  on  char.ncterist ic  (2). 

Finally,  it  was  observed  chat  changes  in  and  0^^^  keeping  the  sum 
equal  to  2000  cal/g,  resulted  in  no  perceptible  change  in  ignition  time. 
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CASE  C 

Case  C Is  similar  to  Case  B in  having  two  rate  processes:  (1)  an 

endotheriDlc  surface  pyrolysis  supplying  fuel  for  the  (2)  exothennic  gas 
p ase  reaction  with  ambient  oxidizer.  The  same  number  of  parameters  is 
required  for  a complete  description.  Qualitatively,  the  surface  pyro- 
lytic process  limits  the  surface  temperature,  but  here,  the  similarity 
with  Case  B ends.  Even  though  the  surface  temperature  is  restricted, 
two  extremes  of  ignition  characteristics  may  be  noted.  First,  a fast 
gas  phase  reaction  occurs  almost  immediately  after  pyrolysis  commences, 
and  second,  a slow  gas  phase  reaction  develops  into  ignition  only  after 
a time  delay  for  diffusional  mixing  and  accompanying  reactive  heating. 
Both  aspects  are  discussed  by  Klndelan  and  Wllliams^^  and  will  not  be 
elaborated  upon  hero.  In  Reference  15,  some  computer  results  are  pre- 
sented for  the  Case  C model.  An  interesting  conclusion  is  that  for  a 
given  pyrolysis  activation  energy  there  is  an  upper  and  a lower  reaction 
activation  energy  outside  of  which  ignition  does  not  occur.  Likewise, 
a given  reaction  energy  has  upper  and  lover  limiting  values  of  pyrolysis 
activation  energy.  In  the  current  study  the  only  effect  noted  (qualita- 
tively) is  that  for  given  activation  energies  of  the  two  reactions,  there 
appear  to  be  limits  on  the  relative  values  of  pre-exponential  factors 
for  which  ignition  occurs.  The  ignition  criterion  was  taken  as  the 
appearance  of  a temperature  maximum  (spatial)  in  the  gas  phase. 

CASE  D 

Case  D,  with  three  chemical  reactions,  i.*  the  most  resistant  to  a 
comprehensive  parameter  study.  It  also  may  represent  the  most  practical 
model,  that  of  gas  phase  ignition  of  two  products  of  surface  pyrolysis. 
Even  a breakdown  and  study  of  the  various  subcases  according  to  the  rela- 
tive values  of  the  three  activation  energies  would  present  a formidaole 
task.  In  the  special  subcase  where  equal  pyrolysis  energies  prevail, 
ignition  certainly  occurs  near  the  surface  at  the  onset  of  pyrolysis  for 
a fast  gas  phase  reaction.  A slow  gas  phase  reaction  would  allow  the 
pyrolyzed  reactants,  already  mixed  by  simultaneity  of  pyrolysis,  to  pro- 
ceed some  distance  and  time  into  the  gas  phase  before  reacting  in  a 
manner  approaching  that  of  a bulk  explosion.  If  one  pyrolysis  product 
has  a much  lower  activation  energy  than  the  »th  r,  it  will  fill  the  gas 
phase,  simultaneously  restricting  the  surface  temperature  so  that  the 
second  Ingredient  will  never  pyrolyze.  There  will  be  no  ignition  in 
this  event-  Just  how  close  together  the  two  pyrolys-.s  activation  ener- 
gies can  be  and  still  permit  ignition  when  both  pyrolyses  are  endothermic 
has  not  been  determined.  The  jnly  systematic  variation  of  parameters  in 
this  study  has  been  the  gas  phase  activation  energy.  Reference  values 
arc  shovTii  in  Table  Including  equal  activation  energies  of  20  kcal/mole 
for  .ill  three  reactions.  Ignition  is  taken  as  the  development  of  a 


Kindel.in,  M. , and  F.  A.  Williams.  ’’Radiant  Ignition  of  a Combus- 
tible Solid  with  Gas-Phase  Exothernicity. ” (To  be  published.) 
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teaperaturc  maxlDun  in  the  gas  phase.  With  a gas  phase  activation  energy 
of  15  kcal/molc.  ignition  occurs  in  10.5  Billiseconds  at  65  laicromcters 
from  the  surface.  This  maximuin  develops  into  a steady  state  temperature 
distribution  after  1.2  additional  milliseconds,  with  the  maximum  moving 
to  33  micrometers  from  the  surface.  At  20  kcal/mole.  a runaway  gas  phase 
teupcrature  occurs  in  31.6  milliseconds  at  70  micrometers  from  the  sur- 
face. At  25  kca'i/mole.  no  ignition  occurred  before  termination  of  the 
run.  steady  state  pyrolysis  having  been  established  at  60  milliseconds. 


TABLE  4.  Reference  Values  for  Case  D. 


Parameters 

Reference  value 

Pyrolysis  reactions: 
Activation  energy 
Heat  of  pyrolysis 
Pre-exponential  factor 

20  kcal/mole 
-150  cal/g 
3 X 10’ 

Gas  phase  reactions: 
Activation  energy 
Heat  of  reaction 
Pre-exponential  factor 

20  kcal/mole 
2300  cal/g^ 
3 X 10“ 

External  heat  flux 

10  cal/cm*  sec 

SUMMARY  AND  DISCUSSION 


Four  simplifications  of  the  unified  Ignition  model  have  been  consid- 
ered for  numerical  study,  two  each  involving  exothermic  chemical  reaction 
at  the  surface  and  in  the  gas  phase.  The  complexity  of  the  particular 
case  has  determined  the  degree  of  completeness  of  the  results. 

Case  A has  been  treated  by  an  extension  of  a previous  analysis  to 
include  gas  phase  convection  and  surface  regression,  both  of  which  were 
foMnd  to  be  important  in  restricting  the  transition  to  steady  state.  A 
unique  feature  of  Case  A is  the  convective  blocking  of  oxidizer  from 
reaching  the  surface  so  that  other  processes  must  augment  a direct  re- 
action between  solid  and  gas  phases  to  achieve  Ignition. 

The  sensitivity  study  on  Case  B is  restricted  to  a single  reference 
point  but  is  complete  otherwise.  While  interactions  may  modify  the 
trends  for  other  reference  points,  activation  energy  is  shown  to  be  the 
most  important  single  parameter  in  the  model.  However,  when  the  kinetics 
is  varied  so  as  to  maintain  constant  reaction  rate,  the  effect  of  acti.a- 
tlon  energy  vanishes.  Variation  of  pressure  and  transport  properties 
has  little  effect  on  ignition  time  but  restricts  the  ignition  regime. 

The  relation  between  ignition  time  and  external  heat  flux  appears  simi- 
lar to  relations  obtained  for  solid  phase  bulk  reactions  and  simple  sur- 
face exothermic  reactions  so  that  such  relations  cannot  be  used  as  a 
diagnostic  tool  for  the  purpose  of  selecting  the  appropriate  model. 
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When  ve  consider  gas  phase  exothermic  processes « We  find  that  under 
certain  conditions,  namely,  those  leading  to  rapid  gas  phase  reaction, 
the  effect  of  external  stimulus  level  on  ignition  time  is  again  similar 
to  results  obtained  for  surface  exotherms*  A unifying  statement  can 
thus  be  made  on  the  basis  of  the  sparse  studies  now  available:  The 

nature  and  site  of  an  exothermic  process  cannot  be  distinguished  by 
q versus  t relationships  if  that  process  occurs  near  the  surface.  The 
effect  of  varying  the  concentration  seems  to  offer  no  solution  because 
of  the  small  effect  of  concentration  versus  temperature  on  reaction  rate. 
Concentration,  as  veil  as  transport  effects,  do  tend  to  impose  limita- 
tions on  ignition  regimes  so  that  calculations  and  experiments  involving 
interrupted  heating  may  offer  a usable  diagnostic  tool. 

For  a slow  gas  phase  reaction,  ignition  times  are  significantly  in- 
creased above  the  thermal  ignition  time  by  dlffusional  effects.  Gas 
phase  reactions  in  this  category  should  be  distinguishable  experimentally 
from  any  surface  reaction  by  the  use  of  high-speed  photographic  techniques. 

Additional  analytical  work  i.  needed  to  settle  the  unresolved  ques- 
tions remaining  from  the  limited  coverage  at  this  time.  All  ansvers  will 
not  be  obtained  by  the  use  of  numerical  techniques  because  of  the  great 
number  of  parameters  involved.  Rather  more  rewarding  should  be  the  de- 
velopment of  ignition  criteria  based  on  reasonable  physical  arguments 
and  mathematical  approximations  and  tested  by  spot  cbecKipg,  using  a 
numerical  method  such  as  the  unified  ignition  theory. 
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NOKENXLATURE 

Theroal  permittivity,  p//c 

Stoichiometric  constants 

Specific  heat 

Mass  diffusivity 

Activation  energy 

Thcmal  enthalpy 

Lewis  number,  pcO/A 

Surface  mass  flux 

Heats  of  reaction  or  pyrolysis 

External  heat  flux 

Universal  ga-;  conrtant 

Linear  burning  rate 

Temperature 

Time 

Mass  fraction 

Diffusion  parameter  (EU]uatlon  23) 
Activation  energy  parameter  (Equation  24) 
Solid  and  gas  phase  responsivitlcs,  ^ 

r,  ^ 

External  heating  parameter  (Equation  25). 

Reduced  mass  fraction 

Reduced  temperature,  (T  - r^)/rQ 

Thermal  conductivity 

Molecular  weight 

Pre-exponent ial  factor 

Diii>ens  ion  less  length 

Density 
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I Dimensionless  time  (different  for  Case  A and  the  ocher  cases) 

^ Hwarth  transformed  length  (/pdx/PQ) 

Subscripts 

0 Initial  value 

1 Solid  phase 

2 Gas  phase 

s Surface 

ss  Steady  state 

couv  With  convection 
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